A Transfer Matrix Study of the Staggered BCSOS Model 
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The phase diagram of the staggered six vertex, or body centered solid on solid model, is inves- 
tigated by transfer matrix and finite size scaling techniques. The phase diagram contains a critical 
region, bounded by a Kosterlitz-Thouless line, and a second order line describing a deconstruction 
transition. In part of the phase diagram the deconstruction line and the Kosterlitz-Thouless line 
approach each other without merging, while the deconstruction changes its critical behaviour from 
Ising-like to a difi'erent universality class. Our model has the same type of symmetries as some 
other two-dimensional models, such as the fully frustrated XY model, and may be important for 
understanding their phase behaviour. The thermal behaviour for weak staggering is intricate. It 
may be relevant for the description of surfaces of ionic crystals of CsCl structure. 



PACS numbers: 68.35.Rh, 64.60. Cn, 64.60.Fr, 6f 
I. INTRODUCTION 

Six vertex models were introduced by Slater to de- 
scribe ferroelectricity in two dimensional networks. Plac- 
ing arrows on the bonds of a square lattice one can define 
the sixteen possible arrangements of arrows pointing to- 
wards and away from a lattice point as vertices. In six 
vertex models only those six vertex configurations are 
kept (see Fig. |^) that satisfy the ice rule, i.e. they have 
two arrows pointing in and two pointing out at each ver- 
tex. Assigning energies ei , . . . , ee to these vertices one 
obtains a class of exactly solved models . 

Six vertex models can also be interpreted as surface 
models, by mapping them to the so-called body centered 
solid on solid (BCSOS) models Q, defined as limiting 
cases of a lattice gas, or Ising model, on a body cen- 
tered cubic lattice. Therefore the phase structure of the 
six vertex model as function of its vertex weights can be 
translated directly to the surface phase structure of the 
corresponding BCSOS model. The mapping turned out 
to be important in understanding the properties of the 
roughening transition ^ . Using the exact solution of the 
six vertex model it was found that roughening is a tran- 
sition of infinite order of Kosterlitz-Thouless (KT) type, 
confirming previous renormalisation group results . 

Experimental situations often are too complex to allow 
even a qualitative description by the exactly solved BC- 
SOS models. Various extensions of the standard six ver- 
tex model have been proposed to deal with these cases. 
Two main classes may be identified: one where inter- 
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actions between vertices are added, and another one in 
which the vertex lattice is split into two sublattices with 
different vertex energies. These modifications, however, 
lead to models which, apart from some particular cases 
0,11 , lose the property of being exactly solvable. Other 
techniques (e.g. numerical ones) have to be adopted. 
Models in the first class have been proposed to account 
for further neighbour interactions between surface atoms, 
which may change the symmetry of the ground state and 
give rise to phase transitions other than the roughening 
transition. Vertex interactions were introduced to repro- 
duce the (2 X 1) reconstruction of the (110) face of fee 
noble metals like Au and Pt ^ . This led to investigations 
on equilibrium phase transitions on these surfaces as well 
as on surfaces of lighter metals like Ag, Rh, etc... A 
model of the (100) surface of an fee crystal exhibiting a 
(2 X 2) reconstructed ground state has recently ex- 
tended the list. The second class of models, with vertex 
weights alternating on the two sublattices, are known as 
staggered six vertex models. A staggering only involving 
the weights of vertices 5 and 6 corresponds to the impo- 
sition of a "staggered field" , i.e. a field coupled to the 
arrow directions that changes sign between neighbouring 
arrows. This gives rise to an inverse roughening transi- 
tion in part of the phase diagram Alternating the 
values for the energies of the vertices 1, 2 and 3, 4 on the 
two sublattices leads to a model known as "the staggered 
six vertex model" (or staggered BCSOS model) in the lit- 
erature. In a large part of its parameter space it can be 
mapped onto the Ashkin- Teller model |jl^. Using this 
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transformation Knops investigated the phase structure 
of the staggered BCSOS model in part of its phase di- 
agram by renormahsation group methods [l^ , but until 
recently a large region of the phase diagram has remained 
unexplored. 

In this paper we present a complete account of our 
investigations, over the full range of parameters, of the 
staggered BCSOS model. A preHminary description has 
been given already in Ref. here we present further 
details as well as a number of new results. 

In the unexplored region of the phase diagram the 
model has a ground state which is twofold degenerate, 
therefore it has a symmetry of Ising type. The twofold 
degeneracy is lost at a second order transition line which 
approaches another line of KT roughening transitions. 
The interplay between the two is particularly interesting, 
especially since a similar interplay between a KT and a 
second order transition has been found for several differ- 
ent models, among which other models for reconstructed 
surfaces [^^^, but also the fully frustrated XY model 
0-|l9| and coupled XY - Ising models They 
have received a great deal of attention in recent years 
and till now their critical behaviour is not fully under- 
stood. The strong interplay between Ising and KT de- 
grees of freedom may lead to several possible scenarios 
where, in a certain region of the phase diagram, either 
the two transitions occur close to each other but remain 
separate, or they merge into a single phase transition, 
which may perhaps belong to a new universality class. 

Apart from these more theoretical aspects the model 
is likely to be relevant for the study of the equilibrium 
properties of a certain class of crystal surfaces, e.g. the 
(001) surface of ionic crystals of the CsCl structure. This 
too will be discussed in some detail. 

The paper is structured as follows. In Section II we 
give a description of the model. In Section III we present 
its full phase diagram. In Section IV we review the tech- 
niques employed in our studies, i.e. the transfer matrix 
method and finite size scaling, and discuss the correlation 
functions and free energies we calculated to derive our re- 
sults. In Section V we discuss the critical exponents of 
the model and some possible scenarios for the changes in 
the critical behaviour along the deconstruction line. In 
Section VI we conclude with a brief discussion of related 
models. 

II. THE STAGGERED SIX VERTEX MODEL 

The partition function of the six vertex model is given 

by 

Z = ^ e-^SLi'^«(C)e» (1) 
{C} 

where the sum runs over the set of all allowed vertex 
configurations {C} and ni{C) denotes the number of ver- 



tices of type i in the configuration C {(3 = l/ksT, with 
ks Boltzmann's constant and T the temperature). The 
model has been solved exactly for any choice of values 
of the energies (z = 1 . . . 6). A relatively simple choice 
of the vertex energies is given by ei = e2 = £3 = £4 = e 
and £5 = eg = which defines, for e > 0, the so-called 
F model. The ground state is twofold degenerate and is 
composed of vertices 5 and 6 arranged alternatingly in 
a chess board configuration. The low temperature phase 
is usually called an "antiferroelectric" phase, since along 
both horizontal and vertical rows the arrows predomi- 
nantly alternate in direction. From the exact solution 
it is known that this system undergoes an infinite or- 
der phase transition to a disordered paraelectric state at 
/3e = In 2. 

As already pointed out in the introduction, the six ver- 
tex models are isomorphic to a class of solid on solid 
(SOS) models called BCSOS models |]. Microscopic 
configurations of an SOS model are given in terms of 
discrete heights hi of surface atoms with respect to a 
reference plane. All lattice sites up to these heights are 
occupied and all sites above them are empty. In the BC- 
SOS model the height variables are placed on the dual 
lattice of the six vertex lattice. This is subdivided into 
an even and an odd sublattice, which are intertwined in 
a chess board pattern and on which the surface heights 
assume even respectively odd values only. The even sites 
will be referred to as black (B) sites and the odd ones 
as white (W) sites. In addition the height differences be- 
tween neighbouring sites are restricted to the values ±1. 
The mapping of six vertex configurations to correspond- 
ing configurations of a BCSOS model is very simple. The 
height differences between neighbouring sites are put in 
a one-to-one correspondence with the arrow directions 
in the six vertex configuration. The convention is that 
the higher of the two surface sites is at the right side of 
the arrow. Given a configuration of vertices, the config- 
uration of heights is fixed uniquely once the height of a 
reference atom has been fixed (see Fig. |l|). 

The vertex energies can be reinterpreted in terms of 
bond energies between the atoms. When periodic bound- 
ary conditions are applied along the (say) horizontal di- 
rection of the vertex lattice the number of vertices 5 and 
6 per row is equal, therefore with no loss of generality one 
can always choose £5 = eg = 0, fixing the point of zero 
energy. The vertices 5 and 6 describe local configura- 
tions in which the height variables on either diagonal are 
equal (see Fig. |l|). Vertices 1, 2, 3 and 4 correspond to 
configurations where the height variables along either of 
the two diagonals are different, therefore ei, £2, £3 and £4 
can be viewed as energies needed to break a next nearest 
neighbour bond and produce a height difference of two 
vertical lattice units between neighbouring sites of equal 
colour. 
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In the ordinary BCSOS model the distinction between 
B and W atoms has been introduced only for convenience 
of description, but the two sublattices are equivalent and 
are treated exactly on the same footing. Knops [Q ex- 
tended the model to a two component system where the 
B and W atoms are physically different. While energy 



zero is still attributed to all vertices 5 and 6, Knops as- 
signed two different energies, e and e', to broken bonds 
between W-W and B-B atoms respectively. In terms of 
the six vertex representation also the vertex lattice is di- 
vided into two alternating sublattices I and II on which 
the vertices assume different energies as follows: 



on sublattice I: ei = £2 = e; £3 = £4 = e'; £5 = ^6 = 
on sublattice II: £i — ^2 — e'; £3 = £4 = e; £5 = ^6 = 



(2) 



This choice defines the staggered six vertex model. In the BCSOS representation the model is described by the 
hamiltonian: 



W iW I 

"■3 I 



{ki) 
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(3) 



subject to the constraint that the height difference be- 
tween neighbouring B and W sites is ±1. The first sum in 
runs over all pairs of neighbouring W sites on the sur- 
face and the second sum over the corresponding B pairs. 
Throughout this article we will also use the parameter 5, 
defined by the relation e' — e + 26. As mentioned already 
in the Introduction the model defined here will be re- 
ferred to as the staggered six vertex (or BCSOS) model. 
Obviously when e = e' {5 — Q) one recovers the usual F 
model. 

For negative values of the vertex energies e and e', the 
system may model ionic crystals of bcc structure as for 
instance CsCl . The constraint of minimal height dif- 
ference between neighbouring surface sites reflects the ef- 
fects of the strong attraction between oppositely charged 
ions, while neighbouring pairs of the same colour, hav- 
ing equal charges, repel each other. It is further assumed 
that on top of the Coulombic repulsion other interactions, 
as for instance spin exchange, generate a slight difference 
in the energies for broken bonds between B-B and W-W 
pairs (e ^ e'). In the staggered BCSOS model the inter- 
action range is limited to next nearest neighbours and to 
have a more realistic representation of ionic crystals one 
needs to extend the interactions to further neighbours. 
Yet we expect the phase structure described here for the 
staggered BCSOS model may be encountered in real ionic 
crystals. 



III. THE PHASE DIAGRAM 

We have investigated the phase diagram of our model 
by means of transfer matrix and finite size scaling tech- 



niques, which will be the subject of Section IV. Here we 
present the main results. Since the model shows a triv- 
ial symmetry upon exchange of e and e', corresponding 
to the replacement (J, e) < — > (— (5, e + 25), we can re- 
strict ourselves to the region 5 > 0. The phase diagram 
naturally divides into three sectors of globally different 
behaviour, though smoothly connected to each other. 
These are described in the three following subsections. 



A. The range e > 0; e' > 



For positive values of e and e' Knops ||T^ investigated 
the phase diagram through a mapping onto the Ashkin- 
Teller model the phase diagram of the latter had 
been obtained before by renormalisation group methods 
p^ . On the /3e axis the (001) surface in the correspond- 
ing BCSOS model is in a flat phase for (3e > In 2, whereas 
the interval /3e < In 2 represents the temperature region 
in which the surface is rough. The infinite order transi- 
tion occurring at the KT point /3e = ln2, /3<5 = corre- 
sponds to the roughening transition of this surface. 

Roughening is a phase transition which can be char- 
acterized by the vanishing of the free energy of a step, 
separating two surface regions of different average height. 
The roughening transition results into a proliferation of 
steps leading to a delocalisation of the surface position 
and to a logarithmic divergence of the mean square height 
difference at large distances: 



G{R^j) = {{h, - hjf) ~ 2alrK{T) \nR,j for i?y 00 and T > 



(4) 
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with av the vertical lattice spacing, K{T) a tempera- 
ture dependent prefactor, Rij the distance between the 
lattice sites i and j and Tr the roughening temperature. 
Below Tn, G{R) saturates for large i? at a temperature 
dependent constant value. Renormalisation group calcu- 
lations 1^ show that at Tr the prefactor assumes the 
universal value 



(5) 



In the particular case of the exactly solved F model, K{T) 
is known for every temperature above Tr [Em , that is 



K{T) 



1 



IT arccos A 



(6) 



where A = 1 — In fact Eqs. (^) and ^ are valid 

not only in the high temperature phase of the F model, 
but also for (5 = 0, /3e < 0, which defines the so-called 
inverted F model |^2| : all along the negative (3e axis the 
surface is in a rough state. 

For 5 ^ Q Knops found two critical lines originating 
from the KT point and running into the regions 5 > Q 
and (5 < 0. The lines represent phase transitions of Ising 
type from an ordered low temperature phase to a disor- 
dered flat (DOF) phase, similar to the phase introduced 
by Rommelse and Den Nijs ||2^ . 

The ground state of the model is twofold degenerate. 
At higher temperatures the more weakly bound sublat- 
tice fluctuates above and below the more strongly bound 
sublattice, which remains almost localised at a given 
level. In the limit <5 — *■ (X) the model can be mapped ex- 
actly onto the two dimensional Ising model, which is criti- 
cal at /3e = ln(l + \/2); the strong sublattice is "frozen" to 
height (say) zero, the only freedom left for the heights of 
the other sublattice is to take the values ±1 just below or 
above that of the strong sublattice (see Fig. Accord- 
ing to the renormalisation group results obtained for the 
Ashkin-Teller model the phase transition remains in 
the universality class of the two dimensional Ising model 
all along the critical line down to <5 = 0. 

Starting from the low temperature phase and increas- 
ing the temperature, the system undergoes an Ising tran- 
sition to the DOF phase, while roughening is pushed up 
to infinite temperature |Q . We reinvestigated this part 
of the phase diagram with the transfer matrix methods 



to be described in Section [V and obtained results in full 



agreement with those of Knops. 



B. The range e < 0; e' > 

When e becomes negative the ground state of the sys- 
tem changes drastically (irrespective of the sign of e'). 
Breaking bonds between white atoms now lowers the en- 
ergy, so that at zero temperature one finds the black sub- 
lattice unbroken (provided 5 > Q) while atoms of the 



white sublattice are found alternatingly above and below 
the black sublattice (see Fig. ||). This surface configura- 
tion is commonly referred to as a c(2 x 2) reconstructed 
surface. In the equivalent six vertex representation the 
ground state is formed by columns of vertical arrows run- 
ning alternately all upwards and all downwards, and by 
rows of horizontal arrows running alternately all right 
and all left. Such an arrangement of directed paths is 
known as a Manhattan lattice^ due to its resemblance to 
the one-way street pattern of Manhattan. As the energy 
is invariant under the reversal of all arrows, the ground 
state is twofold degenerate, just as in the case e > 0. 
Indeed in the limit (5 — s- oo the model can be mapped 
exactly onto an antiferromagnetic Ising model leading to 
the value /3e = — ln(l + -\/2) for the critical temperature. 
This constitutes a horizontal asymptote, as in the case 
e > 0, for a second order transition line, whose existence 
can be deduced again from the mapping of the staggered 
six vertex model onto the Ashkin-Teller model |^^. It 
separates a low temperature c(2 x 2) reconstructed phase 
from a high temperature DOF phase, where the recon- 
struction order is lost but the surface is still globally flat: 
it is the same DOF phase found for e > 0; no singulari- 
ties are met in crossing the /35-axis. Our transfer matrix 
calculations confirm the existence of this critical line and 
show it exhibits Ising type critical behaviour throughout 
sector B. 



C. The range e < 0; e' < 

For e' < the mapping of the staggered six vertex 
model to the Ashkin-Teller model leads to negative Boltz- 
mann weights in the latter. It loses its physical relevance 
and cannot be used any more to make predictions on the 
phase behaviour of the staggered six vertex model. In 
spite of this Kohmoto et al. |29| ] have made some conjec- 
tures, which have proven to be correct, on the physical 
situation beyond the "horizon" e + 2(5 = 0. 

Our transfer matrix analysis shows the existence of 
three phases: a low temperature c(2 x 2) reconstructed 
phase and a DOF phase, which are present already in 
sector B, and a rough phase, which is found only in the 
present sector. Two critical lines separate these phases, 
as shown in Fig. ^ the first one is just the continua- 
tion of the second order line beyond the horizon. It still 
separates the c(2 x 2) region from the DOF region and 
asymptotically approaches the axis (35 = Q. We have 
strong indications that, within at least a major part of 
the sector e' < 0, this line does not belong to the Ising 
universality class. We will present the evidence for this 
in Section The other critical line is a line of KT points 
separating the rough region (or critical fan, as predicted 
already by Kohmoto et al. |^^) from the DOF region. 

The point where the KT line meets the vertical axis can 
be determined from the exact solution of the F model as 
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the point where the prefactor of the logarithmic term in 
the mean square height difference is four times as large as 
its universal value assumed at the ordinary roughening 
temperature of the F model 

KiT) = ^K{Tr) = 4 (7) 

from which one obtains, inverting (||): 

/3e = i In (^2 - \/2^ w -0.2674 . (8) 

The factor four in (|^) stems from the fact that for 6 ^ Q 
the roughening transition is driven by steps of a height 
of two vertical lattice units (as in Ref. due to the 

inequivalence between the two atomic sublattices. 

A simple estimate of the roughening transition temper- 
ature based on a random walk approximation (see Ref. 
£1) yields 

This line has been drawn in Fig. |^. Indeed, for large 
and negative [3e it is seen to run very close to the KT 
line, which we could determine with great accuracy by 
the methods described in the next section. 

A most remarkable feature of our phase diagram is the 
apparent merging of the second order and the KT line 
into a single line (see Fig. |^). 

Their horizontal distance d as a function of /?e can be 
well described by a curve of the form d{l3t) — Ce°'^^'^\ 
with a ~ 12 ijl^. This exponential fit suggests that al- 
though the two lines are coming rapidly closer together 
as Pe is decreasing, they do never actually merge. Other 
fits, of the form d{(3e) = C \(3e — /3eo|", which would be 
expected to work in the case of a merging of the lines 
at e = eo, could not be stabilised against changes in the 
fitting range. 

The apparent non-crossing of the two critical lines at 
first looks very surprising. At low temperatures a domain 
wall between two different Ising phases mainly consists 
of diagonal sequences of vertices 5 and 6, as depicted in 
Fig. ph its energy per unit length approximately equals 
— e/v2. On the other hand a step consists mainly of long 
horizontal and vertical chains of overturned arrows and 
has an approximate energy per unit length of 26. To a 
first approximation steps do not couple with the Ising 
order, since the reconstructed phase remains the same 
at both sides of the step (see Fig. ^) . Hence one would 
expect the KT line (characterized by vanishing step free 
energy) and the Ising line (vanishing Ising domain wall 
free energy) to cross near e + 2^/2 6 = 0. We think that 
the actual non-crossing of the two lines can be explained 
as follows. When temperature is raised, more and more 
closed steps will be formed on the surface as one ap- 
proaches the roughening temperature T^j. On these steps 
the direction of the arrows is reversed. In this way the 



Ising order parameter becomes more and more diluted, 
which will, eventually, strongly reduce the free energy of 
a domain wall. If in the end the closed steps become so 
prolific that they cover on average half of the surface, 
without becoming of infinite length, the Ising order dis- 
appears without roughening of the surface. 

For 26 ^ 0.4|e| the thermal behaviour implied by our 
phase diagram is quite intricate and remarkable. At low 
temperatures the surface is in a c(2 x 2) reconstructed fiat 
phase, then on raising the temperature there is a second 
order transition to a DOF phase, rapidly followed by a 
KT transition to a rough phase. Next there is a reen- 
trant KT transition to the DOF phase. This is an in- 
verted roughening transition similar to the one described 
in Ip^ . Finally, as temperature approaches infinity, the 
system asymptotically approaches a rough phase again. 
Instead for 26 > 0.4|e| the system goes through a single 
phase transition from the ordered to a DOF phase and 
remains flat for all finite temperatures. 



IV. TRANSFER MATRIX AND FINITE SIZE 
SCALING METHODS 

Transfer matrix techniques are frequently used in stud- 
ies of the critical properties of two dimensional systems 
with short range interactions. The construction of the 
transfer matrix (TM) follows a standard procedure and 
the interested reader is referred to the existing literature 
@ for details. 

We use two different transfer matrices, one oriented 
parallel to the axes of the vertex lattice and another one 
tilted of 45° with respect to these axes. We refer to the 
former as vertical TM and to the latter as diagonal TM 
(Fig. 1^ shows a configuration of the diagonal TM). We 
consider a lattice of width N and height M, with pe- 
riodic boundary conditions in the horizontal direction. 
For the vertical TM the subdivision of the lattice into a 
white and a black sublattice, combined with the periodic 
boundary conditions, restricts N to even values. For the 
diagonal TM the horizontal and vertical axis are chosen 
along the diagonals of the vertex lattice and iV can be 
odd as well as even. The element Tij of the matrix is 
defined as the Boltzmann weight of a row of TV vertices 
generated by arrow configurations labeled by the indices 
i and j. One has Ty = if this row of vertices does not 
satisfy the ice rule. For the vertical TM, if i and j are 
identical there are in fact two possible configurations of 
rows of vertices: in this case the transfer matrix simply 
sums their Boltzmann weights. 

There are 2^ different arrow configurations for the ver- 
tical TM, whereas for the diagonal TM this number is 
2^^. The largest values of N we could treat numerically 
were N = 22 for the vertical and iV = 12 for the diago- 
nal TM. Actually, due to the rotation of the lattice over 
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45° the latter should be compared to 12\/2 « 17 for the 
vertical TM. 

In the limit Af ^ oo the partition function per row 
becomes: 



M 



lim {ZNy.Mf'^' = MN) 



(10) 



with Ao(A^) the largest eigenvalue of T, from which the 
free energy per row follows as 



/3/(7V) =-lnAo(A^). 



(11) 



To each state i we associate a polarisation Pi — 
Ni^ — Nil, with Nil ^n-'i -^il total numbers of up 
and down | |30| arrows in the state i. By virtue of the pe- 
riodic boundary conditions in the horizontal direction the 
transfer matrix can be reduced to blocks of fixed polarisa- 
tion, since Ty = if ^ Pj (see, for instance, Ref. 
The so-called central block is the one corresponding to 
zero polarisation and describes a flat surface. The blocks 
with polarisation Pi = ±2 (subcentral blocks) describe a 
surface with a step. 

The difference between the free energies of a surface 
with a step and of a flat surface gives the step free en- 
ergy, which, per unit of length, on an A'^ x oo strip can 
be expressed as: 



Pfs{N) - - (InAi(A^) - InAo(iV)) 



(12) 



where Ao(A^) and Ai(iV) are the largest eigenvalue of 
the central and the subcentral block respectively 
The study of this quantity will allow us to determine the 
roughening temperature. 

The deconstruction transition can be studied by con- 
sidering two correlation lengths, which are both defined 
within the central block. We define the inverse correla- 
tion length as: 



G'W --(lnA2(iV)-lnAo(iV)) 



(13) 



where A2(A^) is the second largest eigenvalue of the cen- 
tral block. The other correlation length can be calcu- 
lated from the diagonal TM as the inverse of the domain 
wall free energy per unit length f^^{N), where fw{N) 
is given by 

/(Ar)- ^'"+^';^'"-" ,14) 

with N odd. Indeed in the diagonal TM an iV x cx3 strip, 
with N odd, is partially frustrated since it cannot accom- 
modate the Manhattan ground state without creating a 
domain wall (see Fig. f^{N) can be interpreted as 
the correlation length connected to the correlation func- 
tion between two disorder variables . 

Conformal invariance p3[ | predicts that, at a critical 
point, the correlation lengths scale as N , so the decon- 
struction transition can be located at the crossing point 



of the curves representing the scaled quantities N/S^d 
and Nj3fw as functions of the temperature for differ- 
ent sizes. In reality, as shown in Figs, ^(a) and (b), 
no perfect crossing is found. Instead, pairs of curves 
obtained for sizes N and N + 2 intersect each other 
in a sequence of points, {PSd{N), PeoiN)), respectively 
{P6w{N), (3ewiN)), which converges to the infinite sys- 
tem critical point {/3Sd, P^d), respectively {P6w , P^w)- 
An extrapolation procedure requiring several iterations 
p5[ is then used to estimate pSn or (3Sw- Of course 
the two independent estimates of the critical point have 
to coincide, which provides a good check on the internal 
consistency and accuracy of our procedures. 

To locate the roughening temperature one has to em- 
ploy a different method. The scaling fs{N) ^ 1/N holds 
not only at the KT transition but also inside the rough 
region, where the surface is in a critical state. There the 
curves Nfs{N), plotted as functions of temperature for 
different values of N, coalesce in the limit N oo and the 
point where they detach from each other can be identified 
as the KT point (see Fig. |9|). For an accurate location 
of Tfj one has to use the universal properties of the KT 
transition which give rise to the scaling prediction p2 36 



NPfs{N) = ^ + ^ 



1 



BlnN 



(15) 



which holds exactly at T = Tr, with A and B nonuni- 
versal constants. The constant 7r/4 is characteristic for 
steps with a height of two vertical lattice spacings. The 
free energy of such a step corresponds to the line tension 
between a vortex-antivortex pair with vorticity 2 in the 
dual representation . The KT transition temperature 
is determined by requiring that a three-point fit of the 
form Nl3fsiN) = Aa + 1/ {A + BlnN) yields Aq = n/A. 
For the extrapolation we used iterated fits in the spirit 
of Ref. . We performed this procedure along different 
lines across the phase diagram, scanning lines with P5 
fixed, lines with (3e fixed and thermal trajectories. 



V. CRITICAL EXPONENTS AND CENTRAL 
CHARGE 

As we noted in the previous section, the critical line 
separating the flat from the rough region can be well 
characterized as a KT line. As we will see, the critical 
properties of the second order line are less well deter- 
mined, especially in the region e' < 0. We will calculate 
critical exponents and central charge pertaining to the 
deconstruction transition using flnite-size scaling meth- 
ods. The two exponents a and v are related to the be- 
haviour of the singular part of the surface free energy 
fsing ^ t^^" and of the domain wall free energy fw ~ 
^ (where t = {T-Td)/Td, Td the deconstruction tem- 
perature). They satisfy the finite size scaling predictions 
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1 d'fjN) 
N dt^ 



and 



N 



dfwjN) 
dt 



N- 



N- 



respectively, valid at the critical point T 
other critical indices we will calculate are 



(16) 



(17) 



Tn . Two 



X — TT- 1™ 7TTT 

2tt w^oo ^d{N) 



T=Tn 



2tt n 



lim NPfwiN) 



(18) 
(19) 



T=Tn 



which represent the exponent of the spin-spin correlation 
function |^l|] and that of the correlation function between 
disorder variables respectively. The numerical errors 
on the values assumed by these quantities are obtained 
as follows. We first evaluate the error on the determi- 
nation of the critical temperature ATd from the quality 
of the extrapolation to iV ^ oo of our finite size data 
p5[ . Subsequently, we extract the values of the expo- 
nents, again by iterated fits, at three different tempera- 
tures: Td - ATd, Td and Td + ATq. This procedure 
allows to determine the maximum possible variation on 
the values of a/i^, v, x and x', thus assigning them an 
error bar. Notice these errors are typically small if the 
critical temperature is determined accurately enough. 

Finally, from conformal invariance ||3^ , |39|] it follows 
that the leading finite size correction to the free energy 
per site of an infinite system with periodic boundary con- 
ditions, /oo, is determined by the central charge (or con- 
formal anomaly) c as 



f{N) 
N 



foo + 



nc 
67V2 



(20) 



In fact we analyzed the central charge using the finite 
size approximation 



ciN,N + 2) = 



3 N^{N + 2) ( f{N) f{N + 2) 



2tt (N + l) \ N 



N ■ 



(21) 



which converges to c in the limit N oo. 

With the techniques described above we find that the 
deconstruction line for e < belongs no doubt to the 
Ising universality class in sector B of our phase diagram. 
Good convergence with increasing size is obtained for the 
critical exponents as well as for the central charge, the 
values of which are: 



a 







1 



1 



1 



8^2 



(22) 



In the region e' < the situation is less clear. The con- 
vergence of the data with increasing system size is worse, 
the values of some of the critical exponents seem to vary 



along the critical line and the central charge cannot be 
determined with any great accuracy. Yet our results seem 
to clearly rule out the possibility that the critical line re- 
mains in the Ising universality class. We present the re- 
sults for the various exponents and for the central charge 
below and then draw some more general conclusions. 



A. The exponent x 



In part of sector C of the phase diagram we find dif- 
ficulties in convergence for the quantities extracted from 
the correlation length Fig. |l^(a) shows the be- 

haviour of the exponent x obtained from Eq. ([l8|), only 
along part of the deconstruction line. The extrapola- 
tion procedure to infinite size is in fact far from triv- 
ial close to the horizon e + 26 = 0, where we find non- 
monotonic behaviour with increasing size for N/^]j(N) 
and even for the sequence f36o{N). In order to give an 
estimate of the exponent x nonetheless, we looked at the 
quantity x{N) = iV/(2<D(^)), where CoiN) is the cor- 
relation length evaluated now at the intersection points 

Fig. ^(b) shows some plots of x{N) vs. N along the 
deconstruction line. The curves 1, 2 and 3 refer to critical 
points in sector B located on the deconstruction line at 
(36 = 0.88, 136 = 0.60 and l36 = 0.45. They show a good 
convergence to the Ising exponent x — 1/8. The other 
curves, (4-10), refer to the values (36 = 0.37, 0.31, 0.28, 
0.25, 0.23, 0.21, 0.20, 0.19 in sector C of the phase dia- 
gram. As the system size increases the curves (4-7) show 
a reentrant behaviour towards the value x = 1/8. At 
values of (36 (S 0.20 we find monotonic convergence again 
as function of system size, but to values which vary con- 
tinuously as shown in Fig. |l^(a). The behaviour of this 
set of curves suggests that along any thermal scan in sec- 
tor C of the phase diagram the quantity x{N) will show 
an asymptotic decrease after a maximum. The position 
of the maximum gradually shifts to higher values of size 
until it exceeds the largest value accessible to our calcu- 
lations and eventally disappears from sight. As already 
mentioned, no accurate fit can be performed on curves 
4, 5, 6 and 7 of Fig. |l^(a), though a rough estimate pro- 
vides values of x below 1/8. When a fitted value can 
be extracted again (at smaller values of ()6) and drawn 
in Fig. |l^(a), one should thus be cautioned against the 
possibility of missing a maximum and a decreasing part. 
This would provide values of x possibly below 1/8 and 
more in accordance with those of x' given in the follow- 
ing sub-section. However, another difficulty may arise: in 
the vicinity of the roughening transition the correlation 
lenght S,D may also be strongly influenced by steplike ex- 
citations . A better quantity to look at is represented 
by the exponent x' . 
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B. The exponent x' 

The quantity N(3fw{N), converges monotonically as 
function of the system size N all along the deconstruc- 
tion line. For 0.25 ^ PS ^ 0.40 the convergence is slow, 
but it is still possible to give an estimate of the expo- 
nent x' using Eq. (^). However the error bars are fairly 
large. We notice a change in the direction of convergence: 
NPfwiN)/2TT converges to x' from above in sector B of 
the phase diagram but from below in sector C. Around 
the line e' = finite size effects are very small. For 
0.3 ^ /3S ^ 0.4 the exponent is still compatible, within 
error bars, with the Ising value of 1/8, as shown in Fig. 
pi] , but for (36 < 0.3 the exponent shifts towards values 
well below this. 



C. The exponents a and v 

Fig. |l2| shows the exponents a and i/ calculated with 
the diagonal transfer matrix along the deconstruction 
line, with the aid of standard extrapolation methods 
based on the scaling relations (|l^,|l^) @. It is al- 
most impossible to obtain these exponents using the ver- 
tical transfer matrix, due to difficulties in convergence 
with increasing size. These problems are much less se- 
vere with the diagonal transfer matrix (see also the Ap- 
pendix), even though the maximum available system size 
is smaller. The values thus obtained are not compatible 
with Ising exponents when e' < 0. They do satisfy the 
hyperscaling relation 2i/ — 2 — a within error bars. 



Ising value c = 1/2. Due to strong finite size effects, slow 
convergence and nearness of the KT line we cannot give 
a reasonable estimate for its actual value. 

Fig. |l^(c) shows the central charge calculated along the 
thermal trajectory e -I- 106 — 0, starting from the rough 
region (at small P6) towards the reconstructed phase at 
larger P6. According to our numerical results the line 
e + 106 — crosses the roughening and the deconstruc- 
tion line in two points very close in temperature. In the 
infinite system limit the central charge should be 1 in 
the rough region, drop abruptly from 1 to at the KT 
point, remain in the DOF region, assume a non-zero 
value at the single point where the trajectory crosses the 
deconstruction line and remain beyond that. In finite 
systems this behaviour is smeared out, as is the case also 
in the other plots in Fig. |l^. Hence, since the KT point 
and the deconstruction transition are extremely close on 
this trajectory, one expects to see an apparent conver- 
gence of c to the sum of the KT value 1 and that of the 
deconstruction transition. For a deconstruction of Ising 
type this would yield c = 3/2. From conformal invari- 
ance it follows that for unitary models with central 
charge smaller than unity c can only assume the values 



1 
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M{M + 1) 



with M = 3.4, 



(23) 



The Ising value c = 1/2 is the lowest possible value, ob- 
tained with AI = 3. Higher values of M correspond to 
phase transitions in different universality classes. From 
Fig. |l3|(c) it is apparent that c converges to a value larger 
than 3/2, which we estimate around c = 1.7 -1.8. 



D. The central charge 

In general the central charge c vanishes in non-critical 
phases (here the flat reconstructed phase and the DOF 
phase) and assumes finite values at critical points or in- 
side critical regions (like the rough phase). As in the 
determination of the exponents a and c is calculated 
with the diagonal TM, as this leads to better convergence 
and smaller finite size effects than calculations with the 
vertical TM. Fig. p^(a) shows finite size approximations 
of c along vertical lines in the phase diagram based on 
Eq. (|2l]). The left part of Fig. |l^(a) refers to a scan 
with (36 — 0.55, which crosses the deconstruction line 
in a point of sector B, where we find exponents in the 
Ising universality class. In this case the central charge at 
the transition shows good convergence towards the Ising 
value (c = 1/2). The right part of Fig. ^(a) refers to 
a scan which crosses the deconstruction line in a point 
of sector C with (36 — 0.25. Fig. |l^(b) shows two other 
plots of central charges along vertical lines with (36 = 0.22 
(left) and (36 = 0.20 (right). In this part of the phase di- 
agram the central charge increases markably beyond the 



E. Deconstruction of non-Ising type? 

The results presented above strongly suggest that the 
deconstruction transition is not in the Ising universality 
class in, at least, part of the region e' < 0. One cannot 
entirely exclude the possibility that the observed devia- 
tions of critical exponents and central charge from their 
Ising values are due to strong cross-over effects, induced 
by the vicinity of the KT line rather than being a 
genuine feature of the deconstruction transition; in view 
of our numerical results however, we believe this is quite 
unlikely. 

Of course, the next intriguing question is: what, if not 
Ising, is the universality class of this reconstruction line? 
The answer to this question is not easy and our numerical 
results are not conclusive. 

In general the exponents vary along the deconstruction 
line, although some vary less than others. The exponent 
X shows generally worse convergence than the exponent 
x' and extrapolation of the values of x in part of the phase 
digram turned out impossible due to the non-monotonic 
behaviour of the finite size data as function of the system 
size N. 
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The exponent x' varies along the deconstruction hne 
as weh, but it remains roughly constant in a limited re- 
gion around the value of /3J 0.2, with small error bars 
thanks to rapid convergence of the finite size data. At 
smaller values of /3(5 its value increases as well as its error 
bars. This may be due to the vicinity of the KT line or 
to the finite size effects caused by the increasing length 
of straight step segments. In general finite size effects in- 
crease at smaller values of /3(5 (see also the Appendix); in 
this part of the phase diagram the most important excita- 
tions consist of closed loops of reversed arrows which may 
become very elongated as the energy per unit of length 
for a straight segment is proportional to 2(5. One should 
expect that finite size effects are particularly strong when 
the typical size of a loop becomes of the same order of 
magnitude as the width of the strip, TV. Slow conver- 
gence also is present in a region to the left of the line 
e' = 0, as can be seen from the large error bars around 
« 0.3 in Fig. |ll|. This is due to a poor determination 
of the value of the deconstruction temperature Td . 

The exponents a and as calculated from Eqs. 
( p^ , p^ , vary along the deconstruction line in sector C. 
However, the hyperscaling relation 2v = 1 — a\s always 
satisfied within error bars. In general, as shown in Fig. 
|T2| , a. tends to have larger error bars than v. In the region 
/?(5 ~ 0.2, the convergence is rapid in the sense that a two 
parameter fit is sufficient to extract a and v from ( p^ ) 
and (p7|). At smaller values of /3(S one in general needs to 
consider corrections to scaling using a three parameter 
fit. 

Unfortunately our numerical results do not allow an 
inequivocal identification of the critical behaviour of the 
deconstruction transition in the region e' < 0. We notice 
however that the exponent x' remains constant in a re- 
gion around /3(5 k, 0.2, where the error bars are smallest. 
In this region also a and v converge rapidly with increas- 
ing size, compared to other parts of the deconstruction 
line in sector C. One possible candidate for the observed 
exponents in this region could be that of the four state 
Potts model, for which a — v — 2/3, compatible with 
our calculated values of a and v. 

Conformally invariant models are classified according 
to the value of their central charge, which can assume 
only discrete values depending on some integer M, as 
given in Eq. (^3|) . At fixed values of M conformal invari- 
ance predicts also the possible values for the expo- 
nents of correlation functions at the critical point. For 
the four state Potts model, the predicted exponents are of 
the type x, x' = 2p^ / with p and q integers, as pointed 
out in Ref. |^l|. For p — 1 and g = 4 one indeed ob- 
tains the well-known magnetic exponent 1/8, instead for 
p = 1, 5 = 5 one obtains the value 2/25. Both values 
are shown as horizontal dashed lines in Fig. ^ the ex- 
ponent x' — 2/25 seems to fit the measured values of the 
exponent very well for [35 w 0.2. For the two dimensional 
Ising model conformal invariance predicts the exponents 



X = 1/8 (magnetic) and x = \ (thermal) only. Thus a 
measured exponent of value x' w 2/25 is a quite clear 
sign of non-Ising critical behaviour. 

Further, the central charge clearly shifts away from its 
Ising value c = 1/2. For the four state Potts model we 
should expect a central charge equal to 1 {M cxi in Eq. 

)). The central charge markably increases in the re- 
gion e' < 0. However, like for the critical exponents, this 
increase goes smoothly from the Ising value, c = 1/2, to- 
wards higher values. The central charge calculated along 
the line e-t- 10(5 = 0, where the deconstruction and rough- 
ening line are almost coinciding in temperature, extrap- 
olates to c = 1.7 — 1.8, well above the Ising plus KT 
value c = 1/2 + 1 — 3/2. As pointed out above, this is 
another indication of non-Ising behaviour of the decon- 
struction transition, though not quite compatible with 
that of the four state Potts model, which would imply a 
central charge equal toc=l + l = 2. 

Finally, also the possibility of having a line with contin- 
uously varying exponents, as the behaviour of especially 
the exponent x in the sector C suggests (Fig. |l0|(a)), 
should be considered. In this case the central charge 
would equal unity, as in the four state Potts model. 

Anyhow, as discussed in Section V A - B, the conver- 
gence of X is much poorer than that of x' . The slow shift 
of x' away from the Ising value as [35 decreases in sector 
C (Fig. ^l]) is known to be a common feature of finite 
size scaling in the vicinity of points where a change of 
universality class occurs. Moreover, as [35 becomes very 
small, the nearness of the KT line is seen to influence the 
convergence of the exponents of the deconstruction line. 
In conclusion, it seems quite plausible to have in practice 
only a window of [35 values where constant critical expo- 
nents are found. To enlarge this window one would have 
to consider larger system sizes. 

Bastiaansen and Knops [pT| recently studied a six ver- 
tex model with an extended range of interactions. They 
also found a phase diagram with a second order line ap- 
proaching a KT line. The exponents of the second or- 
der line clearly deviate from their Ising values and the 
authors suggested they might be explained as tricritical 
Ising exponents. Applied to the staggered BCSOS model 
this would mean a deconstruction line of Ising type with 
a tricritical point, continuing beyond this point as a first 
order line, which is the phase behaviour of the annealed 
diluted Ising (or Blume-Capel |4|]) model. The expo- 
nents at the tricritical point would be a = 8/9, v = 5/9 
and the central charge c = 7/10. For the exponent x' con- 
formal invariance predicts x' = 3/40. Around [35 k, 0.2 
the extrapolated value of x' would also be compatible 
with this value, but a and v are far away from their 
tricritical values. At smaller values of [35 we do find 
exponents which approach those of the tricritical Ising 
model, but this happens in a region where the values we 
obtain for x' clearly shift away from 3/40 and where in 
general finite size effects are quite strong. These same 
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finite size effects also make it impossible to tell whether 
at sufficiently small f3S the deconstruction line becomes 
first order or not. 

The point along the deconstruction line where the 
change of universality class occurs is not sharply deter- 
mined by our numerical results. Wc do not observe an 
abrupt jump of the exponents at a given point, rather a 
continuous shift. A reasonable candidate for the point 
separating the two regions (i.e. Ising and non-Ising), 
could be the point where the deconstruction line crosses 
the line e' = e + 26 = 0. Crossing this line, we find 
changes in the type of convergence of the exponents x 
and x' (Sections V-A, B), although without an abrupt 
change in their values. We recall that in the surface rep- 
resentation of the model, in one region the coupling con- 
stants between the atoms in the two sublattices are both 
negative (e < 0, e' < 0); in the other (where the de- 
construction transition is of Ising type) one of the two 
coupling constants is positive (e' > 0). In terms of the 
vertex lattice, in the region e' > 0, vertices 5 and 6 are the 
excited vertices with the lower energy above the ground 
state value; at e' < 0, vertices 5 and 6 get the higher 
excitation energy. 

All these considerations suggest that the properties of 
the system may change between the two regions e' > 
and e' < and make it more plausible that the shifts 
in the exponents are not just due to cross-over, but also 
result from a real change of universality class of the de- 
construction transition ]4^ ]. 

VI. DISCUSSION AND CONCLUSION 

In this article we studied the critical properties of the 
staggered BCSOS model. Using transfer matrix tech- 
niques we found two critical lines describing the decon- 
struction and the roughening of the (001) surface of a 
two component bcc crystal. 

The two lines approach each other in part of the phase 
diagram, apparently without merging. According to our 
results the deconstruction line in part of the phase di- 
agram changes its universality class from Ising to non- 
Ising, although further investigations are needed to make 
this point more convincing. On the basis of the expo- 
nents we find, we conclude that a possible universality 
class matching these exponents reasonably well, in the 
region where the best convergence is found, is that of 
the four state Potts model. Another possible scenario is 
that proposed by Bastiaansen and Knops [|ll|. In their 
six vertex model with interactions extended to further 
neighbours it is hard too to distinguish between a single 
critical line and two lines approaching each other, but 
remaining separate. They find critical exponents for the 
deconstruction transition clearly deviating from the Ising 
values and conjecture that the observed criticality could 
be explained as tricritical Ising behaviour. The idea of a 



diluted Ising model is particularly attractive in our case 
where, as we have seen, the deconstruction transition is, 
to all likeliness, the consequence of the dilution of the 
Ising order in the system caused by the formation of a 
large number of closed steps of finite length. Unfortu- 
nately we find little numerical evidence for this scenario. 
Finally, also the possibility of having a line with continu- 
ously varying exponents cannot be completely excluded. 
Other models of reconstructed surfaces have been stud- 
ied by several authors. Den Nijs jl^ introduced a model 
that describes (110) missing-row reconstructed surfaces 
of some fee metals {Au, Pt, . . . ). He found a deconstruc- 
tion and a roughening line merging into a single critical 
line, whith Ising and KT behaviour simply superimposed. 
From his data, as presented in the literature, it is not 
possible to really distinguish between actual merging or 
mere rapid approach of the lines. A clear distinction to 
our model is that in Den Nijs' model the deconstruction 
transition remains of Ising type throughout. Another 
class of models for the same metal surfaces has been 
developed and extensively studied by the Trieste group 
||l^,Q. Again a deconstruction line and a KT roughen- 
ing line are seen to approach each other. Depending on 
the microscopic details of the model, the deconstruction 
line keeps its Ising character either all along, or up to a 
tricritical point where it changes to a first order line. 

As mentioned in the Introduction there are several 
other two dimensional models with KT and Ising degrees 
of freedom. One that has received a lot of attention, 
starting from the beginning of the last decade |Q , is the 
fully frustrated XY model, which describes certain two 
dimensional Josephson junction arrays. The study of its 
critical behaviour has led to several different conjectures 
about its universality class and critical exponents. Sev- 
eral papers [ ^t|Ji^ report non-Ising exponents and it has 
been suggested that the model would belong to a novel 
type of universality class. To our knowledge, whether 
this type of universality class would or would not coin- 
cide with that of some known models has not been es- 
tablished yet. In the most recent study concerning the 
fully frustrated XY model Olsson jl^ presents evidence 
of two separate transitions, a KT and an Ising one where 
the former occurs at somewhat lower temperature than 
the latter: Txt < T/s- This would be in agreement 
with our results since the XY model can be mapped onto 
a solid-on-solid model via a duality transformation 
which maps the low temperature phase of one model onto 
the high temperature phase of the other and vice versa. 
Olsson's work suggests that the non-Ising exponents ob- 
served by other groups are due to the failure of some 
finite-size scaling hypothesis used in previous works. In 
the staggered BCSOS model instead, we find clear evi- 
dence of non-Ising exponents. Unfortunately there exists 
no exact mapping between this model and the fully frus- 
trated XY model, therefore they may well be in different 
universality classes. Yet, we hope that some of the ideas 
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developed in this paper to study the staggered BCSOS 
model, will be generalised to other models so as to reach 
a deeper understanding of their critical properties. 
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APPENDIX: FINITE SIZE EFFECTS 

We show here how the diagonal and the vertical trans- 
fer matrix have different finite size effects in part of the 
phase diagram. Consider first the phase point 6 = 0, 
e —oo. Vertices 5 and 6 are absent and the partition 
function can be calculated easily. Consider lattices of 
size N X AI with cylindrical geometry, that is the N ver- 
tices along a horizontal row are connected to each other 
through periodic boundary conditions. In the vertical 
transfer matrix there are in total 2^+*"'^ configurations, 
and the free energy per site is given by: 



g/3e _|_ g 20S _ ^ the vertical and diagonal trans- 



■/?/ 



N + M 
NM 



In 2 



(Al) 



(the symbol / is used to distinguish the free energy per 
site from /, the free energy per row). In the thermody- 
namic limit N,M oo the free energy per site vanishes. 

However in transfer matrix calculations one takes the 
limit Af — > oo keeping N finite; this gives a free energy 
per site equal to: 



(3f = — ln2. 



(A2) 



With the diagonal transfer matrix the total number of 
configurations available is 2^^, since once the arrows on 
a row are fixed the whole configuration is fixed. Repeat- 
ing the same calculation as done above one finds a free 
energy per site: 



-/3/ = 



(A3) 



independent of the value of N. The conclusion is that 
the free energy shows finite size corrections of the or- 
der 1/A^ in the vertical transfer matrix, while there are 
no finite size effects for the diagonal transfer matrix. In 
Fig. n3 we plot the free energy calculated along the line 



fer matrices. The endpoint e 



-2/35 



1 corresponds to 



the free energy which we calculated in (Al) and 
As can be seen from the figure, there is a wide area to 
the left of this point where the free energies calculated 
from the vertical transfer matrix show large finite size 
effects, while the convergence is much faster for the di- 
agonal transfer matrix. In both cases the convergence 
is faster for e^^^^ < 1/2, that is in the region of the 
phase diagram where only a second order line is present. 
Obviously, for small values of 6 the boundary effects are 
very strong, due to closed loops of reversed arrows wind- 
ing around the cylinder, i.e. the ribbon which consti- 
tutes our system with periodic boundary conditions in 
the horizontal direction (see Fig. |l^). These closed loops 
are more frequent in the vertical transfer matrix, since 
one can reverse the arrows along a horizontal line with 
a cost in energy of 2^A^. Closed loops in the diagonal 
transfer matrix require at least a vertex 5 and a vertex 
6, and therefore occur less frequently. 
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FIG. 1. The six vertices and the corresponding height configurations. 

FIG. 2. Side view of the surface for e > in one of its ground states (a) and in the Ising limit /3(5 ~+ oo at a finite value of 
temperature (b). In the ground state atoms form uninterrupted rows also in the [010] direction (orthogonal to the page). The 
structure of the DOF phase resembles that of (b), with broken bonds between B atoms always less frequent than between W 
atoms (their energies are e', e respectively, with e' > e > 0). In the figure, e and e' denote bonds the breakage of which would 
cost these amounts of energy. 

FIG. 3. Side view of one of the ground states of the model for e < 0. Notice that W atoms alternate in height with respect 
to the B sublattice also in the [010] direction (orthogonal to the page). 



FIG. 4. The phase diagram of the staggered BCSOS model. We show here only the sector C, and part of the sector B in 
the inset. Open circles denote the deconstruction line and open squares the roughening line. The estimate for the roughening 
condition, provided by Eq. (9), is shown as a dashed line. Is is almost indistinguishable from the correct curve (squares) for 
/?£ < -1.3. 



FIG. 5. A domain wall separating two different Ising phases for e < 0, (a) in side view and (b) seen from above in the vertex 
lattice. 
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FIG. 6. Side view of a (double height) step as an excitation of the Manhattan ground state. 



FIG. 7. Part of a ground state configuration of an A'^ x oo system with N odd. Due to partial frustration, the system produces 
a domain wall made of a sequence of vertices 5 and 6 (denoted by circles). 

FIG. 8. The scaled correlation lengths N^^^{N) (a) and N/3fw{N) (b) used to analyze the deconstruction transition by 

means of the vertical respectively the diagonal transfer matrix. The curves for different system sizes intersect in a sequence of 
points which for increasing N extrapolate to the deconstruction transition temperature. The figures refer to a scan along the 
thermal trajectory characterized by e/6 = —6.0. The extrapolated values for (a) and (b) (dashed lines) coincide within error 
bars (not shown). 

FIG. 9. The scaled step free energy NPfs{N) for different system sizes along a vertical scan across the phase diagram 
{(36 = 0.14). Coalescence of curves is an indication that the surface is in a rough state. Both the roughening transition and the 
inverse roughening transition are visible: they can be roughly localized in the regions where the curves approach the value 7r/4. 

FIG. 10. (a) The exponent x calculated along the deconstruction line, in the range of values of (35 where the convergence is 
monotonic. (b) Non-monotonic behaviour of x{N) as a function of N: curves 1, 2 and 3 refer to critical points in sector B of the 

phase diagram (and tend to the value 1/2), while all the other curves (4-10) refer to points in sector C (see text). For (4-10), 
due to the limitation in the maximum system size available (larger than that shown here) it is almost impossible to obtain a 
good extrapolation for iV — > oo. 

FIG. 11. The exponent x' calculated along the deconstruction line. The horizontal dashed lines represent the Ising value 1/8 
and the four state Potts exponent 2/25. Error bars smaller than the symbol size are not shown. 

FIG. 12. The critical exponents a and u calculated along the deconstruction line. The dashed line represents their value in 
the four state Potts model. 

FIG. 13. Finite size approximations c{N, N + 2) of the central charge from Eq. (22) (numbers denote the system size A'^) 
across different phases in the phase diagram calculated along vertical scans at ((a), left) (35 = 0.55, ((a), right) (35 = 0.25, ((b), 
left) (35 = 0.22, ((b), right) (35 = 0.20 and (c) along the thermal trajectory e/5 — --10. 0. The inset in (c) is just an enlargement 
of the same graph emphasizing the convergence of the central charge to values larger than 1.5. 

FIG. 14. Free energy along the line e^' + e~^^^ = 1 (Eq. (9)) calculated (a) with the vertical transfer matrix and (b) with 
the diagonal transfer matrix for different system sizes. Notice the different scaling behaviour especially at the point e~^'^^ = 1: 
in (a) the curves for increasing A'^ tend to the infinite system size value (zero) like \n2/N, in (b) they are zero. It is apparent 
that finite size effects are generally smaller in (b). 

FIG. 15. When 5 -C e, finite size effects may show up in the form of closed loops of reversed arrows winding around the 

cylinder. In the diagonal transfer matrix (b) however they are less frequent than in the vertical transfer matrix (a) since they 
require at least one couple of vertices 5 and 6 (circles) for closing up the lace. 



13 




3 



4 




9 



(a) 

0.4 I— 




-2P6 



10 




FIG. 15. 



11 



